Strings in charge-transfer Mott insulators: effects of lattice vibrations and the 

Coulomb interaction 
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Applying the canonical transformation with the 1/A per- 
turbation expansion in the nonadiabatic and intermediate 
regime and the discrete generalisation of Pekar's continuous 
nonlinear equation in the extreme adiabatic regime we show 
that there are no strings in narrow-band ionic insulators due 
to the Frohlich electron-phonon interaction alone. The multi- 
polaron system is a homogeneous state in a wide range of 
physically interesting parameters, no matter how strong cor- 
relations are. At the same time the Frohlich interaction al- 
lows the antiferromagnetic interactions and/or a short-range 
electron-phonon interactions to form short strings in doped 
antiferromagnetic insulators if the static dielectric constant is 
large enough. 

The electron-phonon interaction is strong in ionic 
cuprates and manganites as established both experimen- 
tally |)-| and theoretically f|||. The carriers, doped 
into the Mott insulator, are coupled with the antifer- 
romagnetic background as well. The antiferromagnetic 
interactions are thought to give rise to spin and charge 
segregation (stripes) ||,[l0| . There is growing experimen- 
tal evidence [[Tl]-[l3| that stripes occur in slightly doped 
insulators. Their theoretical studies were restricted so 
far to the repulsive strongly correlated models PJlC|, or 
to an extreme adiabatic limit of the electron-phonon in- 
teraction in narrow |l4|,|l5| and wide band [ p~6|JT7| ] po- 
lar semiconductors and polymers. On the other hand 
there is strong evidence that the nonadiabatic electron- 
phonon interaction and small polarons are involed in the 
physics of stripes [p|Jl2]] . Also the role of the long-range 
Coulomb and Frohlich interactions remains to be prop- 
erly addressed. 

In this letter we prove that the Frohlich electron- 
phonon interaction combined with the direct Coulomb 
repulsion does not lead to charge segregation like strings 
in doped narrow-band insulators, both in the nonadia- 
batic and adiabatic regimes. However, this interaction 
significantly reduces the Coulomb repulsion, which might 
allow much weaker antiferromagnetic and/or short-range 
electron-phonon interactions to segregate charges in 
the doped insulators, as suggested by previous studies 

lira. 

To begin with, we consider a generic Hamiltonian, in- 
cluding, respectively, the kinetic energy of carriers, the 
Frohlich electron-phonon interaction, phonon energy, and 
the Coulomb repulsion as 
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with bare hopping integral f(m), and matrix element of 
the electron-phonon interaction 
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Here i = (m, s), j = (n, s') include site m, n and spin 
s, s' quantum numbers, rii — ejej, Cj, d q are the electron 
(hole) and phonon operators, respectively, and N is the 
number of sites. At large distances ( or small q) one finds 
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and 
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The phonon frequency w q and the static and high- 
frequency dielectric constants in = e^ 1 — £q 1 are 
those of the host insulator [% = c = 1). 

One can apply the canonical transformation and 
the 1/A multi-polaron perturbation theory || to inte- 
grate out phonons, 



S = ^2 Ui [ui(q)d q - H.c] . 

q,i 

The result is |Q 

H = e s He- s = ^ %cj Cj - E p ^ <, 
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where 



era = t(m - n)5 SlS > exp ^^[^(q) - u 3 -(q)]<2q - H.c^j 

(7) 

is the renormalised hopping integral depending on the 
phonon variables, E p — ztX is the polaron level shift and 

n)<5 SiS 'C l t c J +^cj q n,;[u l (q)d q + i?.c] i 

^ Vij =F(m-n)-— 2^7(q) w q cos[q • (m - n)] (8) 
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is the net interaction of polarons comprising the long- 
range Coulomb repulsion and the long-range attrac- 
tion due to ionic lattice deformations. Here A = 
^ q 7(q) 2 a; q /27Vz£ is the dimensionless coupling con- 
stant, t is the nearest neighbour hopping integral and 
z is the coordination lattice number. 

The extention of the deformation surrounding 
(Frohlich) polarons is large, so their deformation fields 
overlap at finite density. However, taking into account 
both the long-range attraction of polarons due to the lat- 
tice deformations and the direct Coulomb repulsion, the 
net long-range interaction is repulsive ||. At distances 
larger than the lattice constant, |m — n| > a = 1, this 
interaction is significantly reduced to 



(9) 
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Optical phonons nearly nullify the bare Coulomb repul- 
sion in ionic solids if eo >> 1, which is normally the case 
in oxides. The kinetic energy term in the exact Hamil- 
tonian, Eq.(6) involves multiphonon events generating a 
residual "polar on- phonon interaction j^] . Below we show 
that in the two opposite limits, the nonadiabatic (w q > t) 
and in the extreme adiabatic ( w q — > 0) regimes, there 
is no charge segregation or any other instability of the 
polaronic liquid due to the Frohlich interaction in doped 
insulators, but only Wigner crystallization at very low 
densities. 

First we consider the nonadiabtic and intermediate 
regime. The properties of a single small polaron with 
the Frohlich electron-phonon interaction were discussed 
a long time ago |l^j2^]. Exact Quantum Monte-Carlo 
simulations [^lj showed that the first order 1/A pertur- 
bation theory is numerically accurate for any coupling 
if the phonon frequency is sufficiently large, u> q > t/2. 
The characteristic frequency of phonons strongly coupled 
with carriers is about w q = 75 meV Q in cuprates, so 
cuprates are in this regime. Hence, one can replace the 
hopping operator in Eq.(6) for its phonon average, reduc- 
ing the problem to narrow-band fermions with the weak 
repulsive interaction, Eq.(9). Next order corrections in 
f/A increase the polaron binding energy with little ef- 
fect on the bandwidth [p2[ . Because the net long-range 
repulsion is relatively weak, the relevant dimensionless 
parameter r s — m*e 2 /eo(47rn/3) 1//3 is not very large in 
doped cuprates. Wigner crystallization appears around 
r s ~ 100 or larger, which corresponds to the atomic den- 
sity of polarons n < 1 -6 with eo = 30 and the polaronic 
mass m* = 5m e typical for cuprates and manganites. 
This estimate shows that small polarons in cuprates and 
manganites are in the homogeneous state at physically 
interesting densities. 

In the opposite adiabatic limit one can apply a discrete 
version of the continuos nonlinear equation p3| proposed 
in Ref. 124] for the Holstein (molecular) model of the 



electron-phonon interaction and extended to the case of 
the deformation and Frohlich interactions in Ref. Jl4],|l5| . 
Applying the Hartree approximation for the Coulomb re- 
pulsion, the single-particle wave- function, ip n (the ampli- 
tude of the Wannier state |n)) obeys the following equa- 
tion 

- t(m)[ip n - ^ n+m ] - ec/) n ip n = Ei/j n . (10) 

m=£0 

The potential 4>n,k acting on a fermion k at the site n 
is created by the polarization of the lattice </> l n k and by 
the Coulomb repulsion with the other M — 1 fermions, 

<f>n,k = <fc + C.fc- (11) 

Both potentials satisfy the descrete Poisson equation as 
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and 
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with A0 n = ^2 m {4>n — 0n+m)- Differently from Ref. 
]l5[ we include the Coulomb interaction in Pekar's func- 
tional J |Q , describing the total energy, in a selfconsis- 
tent manner using the Hartree approximation, so that 

n,p,m^0 

£ |Vn,p| 2 A-Vm,,| 2 
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If we assume, following Ref. |L4j that the single-particle 
function of a fermion trapped in a string of the length 
N is a simple exponent, ip n — iV -1 / 2 exp(i/cn) with 
the periodic boundary conditions, then the functional 
J is expressed as J = T + U, where T = —2t(N — 
1) s'm(irM/N)/[N sin(7r/7V)] is the kinetic energy (for an 
odd number M of spinless fermions) , proportional to 
t, and 



U = -—M 2 I N + — M{M - 1)I N , 



(15) 



corresponds to the polarisation and the Coulomb ener- 
gies. Here the integral In is given by 



In = 
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dx I dy 
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It has the following asymptotics, Fig.l, 



In = 



1.31 + In N 

N ' 
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which is also derived analyticially at large N by the use 
of the fact that sin(N x / 2) 2 / (2tt N sin(x/2) 2 ) can be re- 
placed by a S- function. If we split the first (attrac- 
tive) term in Eq.(15) into two parts by replacing M 2 
for M + M[M — 1), then it becomes clear that the 
net interaction between polarons remains repulsive in 
the adiabatic regime as well because k > eoo. Hence, 
there are no strings within the Hartree approximation 
for the Coulomb interaction. Strong correlations do not 
change this conclusion. Indeed, if we take the Coulomb 
energy of spinless one-dimensional fermions comprising 
both Hartree and exchange terms as [E5| 



E c = 



e 2 M(M - 1) 
NeZ, 



0.916 + In M], 



(18) 



the polarisation and Coulomb energy per particle be- 
comes (for large M >> 1) 



U/M 



e 2 M 



[0.916 + In Af — a(1.31 + In TV)], (19) 



where a = 1 — eoo/eo < 1. Minimising this energy with 
respect to the length of the string N we find 



N = M 1/a exp(-0.31 + 0.916/a), 



(20) 



and 



(U/M) min = - — M 1 " 1 /" exp(0.31 - 0.916/a). (21) 

K 

Hence, the potential energy per particle increases with 
the number of particles so that the energy of M well 
separated polarons is lower than the energy of polarons 
trapped in a string no matter correlated or not. The 
opposite conclusion of Ref. ]lq| originates in an incor- 
rect approximation of the integral In oc N 0A5 /N. The 
correct asymptotic result is Zjv = \n(N)/N. 

One can argue |27|] that a finite kinetic energy (t) can 
stabilise a string of a finite length. Unfortunately, this 
is not correct either. We performed exact (numerical) 
calculations of the total energy E(M, N) of M spinless 
fermions in a string of the length N including both kinetic 
and potential energy with the typical values of = 5 
and eo = 30. The local energy minima (per particle) in 
the string of the length 1 < N < 69 containing M < N/2 
particles are presented in the Table. Strings with the 
even fermion numbers carry a finite current and hence the 
local minima are found for odd M. In the extreme wide 
band regime with t as large as 1 eV the global string en- 
ergy minimum is found at M = 3, N = 25 (E = -2.1167 



eV), and at M = 3,N = 13 for t = 0.5 eV (E = -1.2138 
eV). However, this is not the ground state energy in both 
cases. The energy of well separated d > 2-dimensional 
polarons is well below, less than — 2dt per particle (i.e. 
—6 eV in the first case and —3 eV in the second one in the 
three dimensional cubic lattice, and —4 eV and —2 eV, 
respectively, in the two-dimensional square lattice). This 
argument is applied for any values of eo,eoo and t. As 
a result we have proved that strings are impossible with 
the Frohlich interaction alone contrary to the erroneous 
Ref. @|7). 

The Frohlich interaction is, of course, not the only 
electron-phonon interaction in ionic solids. As discussed 
in Ref. |J], any short range electron-phonon interac- 
tion, like, for example, the Jahn- Teller (JT) distortion 
can overcome the residual weak repulsion of Frohlich 
polarons to form small bipolarons. At large distances 
small nonadibatic bipolarons weakly repel each other 
due to the long-range Coulomb interaction, four times 
of that of polarons, Eq.(9). Hence, they form a liq- 
uid state ||, or bipolaronic-polaronic crystal-like struc- 
tures [^U depending on their effective mass and den- 
sity. The fact, that the Frohlich interaction almost nul- 
lifies the Coulomb repulsion in oxides justifies the use 
of the Holstein-Hubbard model |||9| . The ground state 
of the ID Holstein-Hubbard model is a liquid of inter- 
site bipolarons with a significantly reduced mass (com- 
pared with the on-site bipolaron) as shown recently [ p0[ . 
The bound states of three or more polarons arc not 
stable in this model, thus ruling out phase separation. 
However, the situation might be different if the antifer- 
romagnetic [9|]l0|] and JT interaction |3l|] or any short 
(but finite)-range electron-phonon interaction are strong 
enough. Due to long-range nature of the Coulomb re- 
pulsion the length of a string should be finite (see, also 
Ref. [ [l2l]l4ll ). One can readely estimate its length by the 
use of Eq.(8) for any type of the short-range electron- 
phonon interaction. If, for example, we take dispersive 
phonons, ujq = luq + Slo(cos q x + cos q y + cos q z ) with a q- 
independent matrix element 7(q) = 7, we obtain a short- 
range polaron-polaron attraction as 



v att (n - m) = -E att (5u/u)5\ a _ t 



(22) 



where E att — ~/ 2 luq/2. Taking into account the long- 
range repulsion as well, Eq.(9), the potential energy of 
the string with M — N polarons becomes 



U = —N Z I N 
(0 



NEattSu 



(23) 



Minimization of this energy yields the length of the string 
as 



N = cxp I s 2.31 



(24) 



Actually, this expression provides a fair estimate of the 
string length for any kind of attraction (not only gen- 
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erated by phonon dispersion), but also for the antifer- 
romagnetic exchange and/or Jahn- Teller type of interac- 
tions |32). Due to the numerical coefficient in the expo- 
nent in Eq.(24) one can expect only short strings (if any) 
with the realistic values of E att (about a few hundreds 
millivolts), and the static dielectric constant eo < 100. 

We conclude that there are no strings in ionic doped 
insulators with the Frohlich interaction alone. Depend- 
ing on their density and mass polarons remain in a liquid 
state or Wigner crystal. On the other hand the short- 
range electron-phonon and/or antiferromagnetic inter- 
actions might provide a liquid bipolaronic state and/or 
charge segregation (strings of a finite length) because the 
long-range Frohlich interaction significantly reduces the 
Coulomb repulsion in highly polarizable ionic insulators. 
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Figure captions 
Fig.l. The polarisation energy of small Frohlich po- 
larons trapped in a string depends on its length as 
\n(N)/N. 
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